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Abstract 

The three-dimensional potential equation, motivated by representations of quantum me¬ 
chanics, is investigated in four different scenarios: (i) In the usual Euclidean space E 3 
where the potential is singular but invariant under the continuous inhomogeneous or¬ 
thogonal group 10(3). The invariance under the translation subgroup is compared to 
the corresponding unitary transformation in the Schrodinger representation of quantum 
mechanics. This scenario is well known but serves as a reference point for the other 
scenarios, (ii) Next, the discrete potential equation as a partial difference equation in a 
three-dimensional lattice space is studied. In this arena the potential is non-singular but 
invariance under 10(3) is broken. This is the usual picture of lattice theories and nu¬ 
merical approximations, (iii) Next we study the six-dimensional continuous phase space. 
Here a phase space representation of quantum mechanics is utilized. The resulting poten¬ 
tial is singular but possesses invariance under 10(3). (iv) Finally, the potential is derived 
from the discrete phase space representation of quantum mechanics, which is shown to 
be an exact representation of quantum mechanics. The potential function here is both 
non-singular and possesses invariance under 10(3), and this is proved via the unitary 
transformations of quantum mechanics in this representation. 
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1 Introduction 

It is well known that quantum mechanics may be represented in terms of phase space variables 
ifn. This formulation of quantum mechanics has a history dating back to Wigner ||2l and 
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Weyl ||3l. The analog of the “wave function” in this formulation is a probability distribution 
function in phase space known as the Wigner function. It is also known that one may describe 
quantum mechanics in terms of discrete quantities [4J, [51. Finally, the discrete phase space 
representation of quantum mechanics has been established in references |61 and Q. The 
generalization of this representation to relativistic quantum mechanics and field theory has 
been presented in O, |§^1, and lUOll . In ifTOl quantum electrodynamics was formulated in 
the arena of three-dimensional discrete phase space and continuous time. (Such mixing of 
discrete and continuous operations is mathematically rigorous and has a long history. For 
example, see |[TT1l - |[20l {^ The great advantage of this representation is that the S-matrix 
elements turn out to be divergence free. The mathematical reason for the elimination of the 
divergences is the fact that partial difference-differential equations representing free fields 
have non-singular Green’s functions. This discrefe phase space represenfafion is relatively 
new, and has nof been sfudied in any defail in fhe imporfanf physical arenas of non-relafivisfic 
quanfum mechanics and pofenfial fheory, save for specifically fhe harmonic oscillator Q. 
This is fhe primary mofivafion for fhe currenf sfudy. 

In fhis paper we concenfrafe on fhe partial difference equations represenfing fhe pofen- 
fial equafion in fhe fhree-dimensional discrete phase space (which is a subsef of fhe six¬ 
dimensional confinuous phase space ), and sfudy fhe resulfing non-singular pofenfial, which 
replaces fhe usual singular Coulomb pofenfial. Before deriving fhe resulf in discrefe phase 
space we discuss briefly in secfionj^fhe usual ellipfic pofenfial equafion in fhe differentiable 
manifold E 3 . The Green’s functions and fhe potential function Vb( 2 ;^, x^) exhibif well 

known singularifies. Oufside fhe singularify, fhe pofenfial function is invarianf under fhe 
confinuous group XO(3). Alfhough fhis facf is well known, we will utilize a proof of fhis 
invariance from fhe unifary fransformafion of posifion and momenfum operafors, and P^., 
in quanfum mechanics according fo fhe Schrodinger represenfafion. This is done as a refer¬ 
ence for fhe sfudy of quanfum mechanics and pofenfial fheory in fhe ofher represenfafions in 
subsequenf sections. 

In section we touch upon discrete pofenfial fheory II 2 TII in a fhree-dimensional laffice 
space. Physical applications of such a discrefe pofenfial occur in random walks on a laffice 
II 22 II and in finife elecfrical nefworks 1231 . The represenfafion here is also often ufilized in 
numerical approximations. This discrefe pofenfial is non-singular. However, if does nof admit 
invariance under the continuous group XO(3). Also, this difference Laplacian is not capable 
of providing an exact representation of quantum mechanics. 

In section 1^ we deal with the representation of quantum mechanics in the six-dimensional 
continuous phase space manifold |T]. (Originally, Weyl and Wigner suggested such represen¬ 
tations frj [ 23 .) In this section we also formulate the potential equation in continuous phase 
space , again by using the corresponding quantum mechanics operators Q* and The po¬ 
tential here is singular. However, it does admit the invariance under the group XO(3) as a 
subgroup of canonical transformations. 

In section we investigate the potential theory in the new discrete phase space sce¬ 
nario. Firstly, quantum mechanical operators Q* and Pj- are formulated in this representa- 

O 

tion. These are provided by new finite difference operators which we denote as (5®-^ A and 
—iA^. These are capable of providing an exact representation of quantum mechanics, not 
just a discrete approximation of the continuum theory. The corresponding Green’s functions 
are non-singular and the discrete phase space potential, W{'n}, nf, nf) is analogous to the 

'The particular representation in (H-Go] also respects Lorentz symmetry, as was shown in (9). One confusion 
comes from the fact that it is tempting to think that the spacetime is discrete in these theories, but it is not. 
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potential x^) inEs. The potential equation governing the function TU(n^, n^, n^) 

is proved to be invariant under the continuous group XO(3). The proof here is also based 
on the unitary transformation involving the corresponding quantum mechanical operators Q* 
and Pfc. We also study the particular solution VU(o)(n^, n?, rfi), analogous to the particular 
potential of section(corresponding to the “source” point at the origin). It 

is shown that n?, n?) is non-singular and, when analysis is restricted to an axis, 

VU(o) (0) 0) 2n^) can be expressed in terms of gamma functions. 

We present the variational derivation of second order partial difference equations and the 
equations for boundary variations in an appendix for completion |'8l|. 

We choose units such that ^ = 1. A new constant, £ > 0, will also appear in subsequent 
calculations, which is a characteristic length. Except in certain instances where it is required 
for clarity, we set ^ = 1. 


2 The potential function in the Euclidean space E 3 


In this section we briefly review the usual potential function in E3. Most of the analysis in 
this section is well known but this section serves to set the stage for subsequent results for 
comparisons. Before proceeding we establish the notation as follows. Roman indices take 
values from {1, 2, 3} and summation over repeated indices is implied. Hatted quantities 
refer to transformed versions of the corresponding unhatted quantities. The ordered triple 
(X , X , X G denotes a Cartesian coordinate chart for E3. 

The potential function is a solution to the well known equation 


VV(x\ x^) := x2, x3) = P’^djdkV{x\ x^ x^) = 0. (2.1) 

A particular Green’s function for the elliptic partial differential equation ( |2.1| ) is provided by 

exp [ikj (x^ ~ 


G(x^, x^, x^; x^, x^, x^) 


POO POO POO 

J —00 J —00 J — 


(2 

1 


-00 [{hy + {k2y + {ksY] 

1 


dvr [(a;l _ £1)2 + ( 3.2 _ £2)2 + ( 3.3 _ £3)2]1/2 


dkidk2dk‘i 

( 2 . 2 ) 


In electrostatic phenomena, the physical interpretation of ( |2.2| ) is that it represents the Coulomb 
potential at field point ( 3/ ^ tJC ^ tJC y due to a positive unit electric charge located at source point 
(x^, x^, x^) (in Heaviside-Lorentz units). 

By taking the “source point” at the origin, we can write a particular solution of the poten¬ 
tial equation as 


V(o)(a;\ x^ x^) 


G(x^, x^, x^; 0,0,0) 


1 1 

4vr (x^)^ + (x‘^y + {x^y 


for {x^y + (x^)^ + {x^y > 0, (2.3i) 


lim |V(o)(x^, x^, x^)| =0, (2.3ii) 

\x^ \ —> 00 

, lim |V(o)(a:S x^ x^)| ^cx). (2.3iii) 

-> (0, 0, 0) 

The equation ( |2.3iii| ) demonstrates the well known fact that the above potential diverges in the 
limit (x^, x^, x^) —)• (0, 0, 0). 
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The potential equation (2.1 1 may be variationally derived from an appropriate action (the 
Dirichlet integral 1^ . ^3): 


J[V] := dx^dx^dx^ L(y; yi, y2, ^ 3 ) 1 ^ = v{...),b, = a,v{ 


IDs 




L{y,yi,y2,y3)... ■=^ 


yjyk 


I- 


(2.4i) 

(2.4ii) 


Here, D 3 is a closed, simply-connected, proper subset of E3. The varied function is chosen 
to be 

y :=V (x^, x^) -|- eh{x^,x^,x ^), (2.5) 

with e 7 ^ 0 and h being an arbitrary bounded function of class C^{D 3 -, M). The variational 
principle applied to ( |2.4i I using ( |2.4ii 1 implies that 

6 ^^djdkV{x^, x^, x^) = 0 , ( 2 . 6 i) 

and 


IdDs L 


h{x^,x‘^,x^)6^^dkV{x^, x^, x^)nj{x^,x‘^,x^] 


J|... 


d^s = 0. 


( 2 . 6 ii) 


The equation (2.6il is the Euler-Lagrange equation in D3 C M'^, and the equation (2.6iil is 
the boundary term due to the variation on dD^. In the case of Dirichlet boundary condi¬ 
tions /i(x^, x^, x^)g 75 ^ = 0, and (2.6iil is identically satisfied. Moreover, the only prescribed 
Neumann condition n^djV (x^, x^, = 0 is variationally admissible. 

Now we shall briefly deal with the symmetry group admitted by Laplace’s equation (2.1 1 . 
It is the six-parameter inhomogeneous orthogonal group XO(3). This Lie group is explicitly 
provided by the transformation 


x^ =cP -I- , 

(2.7i) 

=^ab , 

(2.7ii) 

aV\ =r\a\ = 5\ , 

(2.7iii) 

x^ =aj (yff — cP) , 

(2.7iv) 


with the six independent parameters d> and comprising the linear transformation. 

The potential function V (x^, x^, x^), and the potential equation (2.1 1 , treated tensorially 
under (|2.7i[-iv), transform as 


V(x^,x^,x^) =V(x^, x^, x^), 


( 2 . 8 i) 


6 ^^didjV(x^,x^,x^) =S^^a^^a'^^d/jdcV(x^, x^, x^) = 6 '^^dbdcV{x^, x^, x^) = 0 . ( 2 . 8 ii) 


The equation (2. 8 ii I illus trates the invariance of the potential equation p.l| ) under the group 

\ and = (5*^ 


X(9(3) expressed in (2.7i iv). If one sets a\ = 5*^ and = (5*^ the transformations comprise 


the three-parameter Abelian subgroup of translations: 

U(x^ + c^,x‘^ + (?,x‘^ + c^) = U(x^, x^, x^), 


(2.9i) 



or, V{x^,x^,x^)=V{x^-c^,x^-c^,x^ 


(2.9ii) 
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It is known that the potential function V{x^, x^) (also called a harmonic function), 

locally admits the Taylor series expansion in a star-shaped domain of 1^ . Therefore, we 
obtain from ( |2.9ii| ). 


V{x^,x^,x^) =U(x^,x^,x^) 


1 ^2 


+E 


{-l)i 


i=i 


E-Et 

ii=\ 

{h+...+ij=j) 






(9x*i • • • c)x*J 


V{x\ x\ x^) 


(2.10i) 


or. 


U(x^,x^,x^) =exp 


-c^dk U(x\x^, x^), 


6°'’^dadbV{x^ ,x^ ^x"^) =exp 


1 ^2 S', 


-c^dk 


6^'>dadbV{x\ x^ x") 


1 ^2 S', 


= 0 . 


(2.10ii) 

(2.10iii) 


Thus, we have an alternate proof of the invariance of the potential equation ( |2.1[ ) under the 
subgroup of translations. This proof has a quantum mechanical aspect to it and will be also 
useful for subsequent analysis of quantum mechanics in the discrete phase space representa¬ 
tion. 

In the Schrddinger representation of quantum mechanics (also known as wave mechan¬ 
ics) the self-adjoint coordinate and momentum operators acting on Hilbert space vectors are 
provided by 


Pkif := - idk'ip{x^,x'^,x^) , 

(2. Hi) 

:=x^ijj{x^,x‘^,x^) , 

(2.11ii) 

[PkQ^ - Q^Pk] -if = - i6\ ipix^.x^.x^). 

(2.11iii) 

Now consider an unusual quantum mechanical equation 

-5^^PjPk^p =0, 

(2.12i) 

or. 

S^^djdktpix^, x‘^, x^) =0. 

(2.12ii) 


This corresponds to a free particle with no kinetic or potential energy. The equation (2.12iii 
implies that both Re('!/)) and Im('!/)) satisfy the potential equation (2.1 1 . The above equations 
can also be understood as the non-relativistic limit of a massless free particle in relativistic 
quantum mechanics via 





^ = 0 , 


(2.13i) 


or, 

— 6^^djdk'4>{x^,x‘^ 3;^, x^, x^) = 0 . (2.13ii) 
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Here c is the speed of light and the fourth coordinate, is ct. In the non-relativistic limit 


oo, the above relativistic equation reduces to (2.12ii i. 


In the Schrddinger picture, the transformations of Hilbert space vectors and position and 
momentum operators, induced by the translation subgroup of XO(3) are specified by 


P —P 

^3 ~-^3 1 




r/j = exp 


-ic^P, 


i’- 


By equations ( |2.1 li| -iii), the transformations above yield 

= exp [—c^dj] x^, x^). 


(2.14i) 


(2.14ii) 


(2.15) 


The equation (2.151 above is the exact replica of ( |2.10ii| ). Hence the equivalence of such 
transformations on the potential V (x^, x^, x^) and on the zero-energy non-relativistic quan¬ 
tum mechanical particle. We repeat here that most of the analysis in this section is well 
known, but it serves as a reference for the analysis of the potential equation and quantum 
mechanics in discrete representations and phase space . We will also comment on advantages 
and disadvantages of these other representations in the appropriate sections. 


3 Discrete space and the potential function 


There is a characteristic length i > 0 here. (Recall that we use units such that i = 1. This is 
common in numerical analysis.) There are also three discrete variables G {..., —1,0,1, 2...} 


Z. We define various partial difference operators as follows: 

The right partial difference operator is defined as 

Aj/(n\n^,n^) :=/(•••, -f 1, •••)-/(•••, ^^ •••) ■ (3-1) 

The left difference operator is defined as 

A'/(n\n2,n^) :=/(••• , •••)-/(•• • - 1, - • •) • (3-2) 

The mean partial difference operator is defined as 

Aj/(n\n^,n^) := ^ [/(••• ,n^ + I, ■■■)- f{-■■ - I, ■■ ■)] ■ (3-3) 

The discrete potential equation is taken to be ll2T1l 

n^) = 0 . (3.4) 


It is usually defined this way as this combination of discrete derivatives yield a self-adjoint 
potential difference equation. As well, it is variationally derivable in this form ll26ll and this 
is the usual Laplacian of discrete numerical analysis. 

The partial difference equation in discrete analysis, Q, is an analogue of the usual 
potential equation (|2.1|) in the continuum. The discrete potential equation (|3.4[) implies that 


U{'n}^v?,n^) =^| \U{n^ + 1 , v?,n^) + U{n^ — 1 , n^,n^)] 

+ [f7(n^,n^ -f l,n^) -|- U{n^,n^ — l,n^)] 

-|- [U{n^,n^,n^ + 1 ) + U{'n}, it ?— 1 )] | . 


(3.5) 
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Thus, the value U (n^, n^) of the discrete potential (or discrete harmonic function) is the 
arithmetic mean of its values at the six neighboring points of (n^, n^, v?). This property is 
analogous to the mean-value theorem for the harmonic function of the con¬ 

tinuous variables ( tC ^ tlj ^ oc studied in the preceding section. A Green’s function for the 
discrete potential equation ( |3.4| ) is furnished by 

^ l-IT nw nw gjjp 


4(27r) 


f —TT J —TT t/ —TT 


Sin 


^ ) -I- Sim 






Ydkidk2dks , (3.6i) 


(3.6ii) 


<5*-^ AiA' G(- • • ) = -<5nlrii<Jn2n2 5n3n3 • 

A discrete potential function U(^o'j{n^,n‘^,n^), analogous to V(o)(3^^) of (2.3il is pro¬ 

vided by 

U(^Q){n^,n‘^,n^) := G(n^, n^, 0,0, 0) 


exp [i/cj 




4(27r) 


f —TT J —TT J —TT 


Sin 


-p sin"' f 


(^) 


■dkidk2dk3. (3.7) 


^j+sin" 

An interesting property of U(Q'^{'n},n?, rfi) is 

C/(o)(0,0,0) [l8 + 12x/2 - lOx/3 - 7x/6] { 0^ AT [(2 - x/3)(x/3 - x/2)] }^ 

(3.8i) 

|[/(o)(0,0,0)| < oo . (3.8ii) 

Here, the function K[- ■ ■] denotes the complete elliptic integral. Thus, the function ?7(o) (n^, n^, 
in (3.71 is everywhere non-singular. For large values of this function possesses the follow¬ 
ing asymptotic behavior 

12 3 1 , 1 r5[(ni)4 + (n2)4 + (n3)4] 

(7(o)(n ,n ,n ) = 


-P 


47r||n|| 327r||n 


lim |tf(0)(«-)| = 0, 


n 


- 3 } + 0 


1 


n 


(3.9i) 

(3.9ii) 


with ||n|| := y^(nf)2^p7r74j2^F(n^. From the expression in braces in (3.9il it can be 
inferred that the potential violates the symmetry group 0(3). Therefore, although this discrete 
potential remedies the singularity, it is not invariant under the continuous rotation group. 
This, of course, is due to the fact that the underlying lattice structure of the theory does not 
possess rotational invariance. The invariance group of the discrete potential equation ( |3.4| ) is 
a crystallographic space group GTII . 


We finish this section with a brief comment that the difference operators in (3.1 l-(3.3l 
are not capable of providing an exact representation of quantum mechanics, although an ap¬ 
proximation in the form of the discrete Schrddinger equation can be provided 1281. The 
discrete Schrddinger equation has been shown to to have applicability in tight-binding mod¬ 
els for crystal solids Il29l and highly localized phases. The new representation of section 
yields non-singular Green’s functions and an exact representation of the algebra of quantum 
mechanics. 
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4 Continuous phase space and the potential function 

This section relies heavily on the representations of quantum mechanics. The usual Schrodinger 
representation was briefly introduced in equations ( |2.11i i-(2.11iiii. Although not common, 
quantum mechanics can also be formulated in phase space |fT]. We provide here a continuous 
phase space representation of quantum mechanics in the following: 

1 


Pj'fp :=^ {pj - idqj) V’ , (4.ii) 

+ ^9pk) V' {q^,q‘^,q^-,pi,P2,P3) , 




ip = -ism'll; Pi,P2,P3) ■ 

The eigenfunctions of the momentum operator are given by: 

'ip{k){q^,q^,q^]Pi,P2,P3) =^(Pi,P2,P3)exp i{V2kj -Pj)q^ 

with 

Pj-^ik) =kj^(k) ■ 


(4.1ii) 
(4.1 hi) 

(4.2i) 

(4.2ii) 


Here, A{pi,p 2 ,P 3 ) / 0 is an arbitrary function. Similarly, the eigenfunctions of the position 
operator are furnished by: 


'iPix){q^,q‘^,q^',pi,P2,P3) = ^(g^g^g^)exp -i{V2x^ - q^)pj 

with B{q^, q^, q^) ^ 0 being another arbitrary function. 

The potential equation, according to (2.12ii and (|4.1i[ii) is provided by 


(4.3) 


-P’^PPki^ =0, 


or. 


2 L 


( 5 ^^ {pj - idqj) {pk - idqk) ^{q^,q^,q^-,Pi,P2,P3) = 0 . 


(4.4i) 


(4.4ii) 


A special solution, which is the analogue in this representation of the potential function in 

exp [i [V2kj - Pj) q^] 


(2.3i iii) is given by 

-.1 3. 


roo roo roo 


^io)iq ,q ,q ',pi,P2,P3) 


' —00 J —00 


1-00 [{kl)^ + (^2)2 + (ks)^] 


dkidk2dk‘i , 

(4.5i) 


or. 






lim |Q(o)(---)| =0, 

\\q\\^oo 


lim 

l|q|Ko+ 


l«(0)('")| 


00. 


for 11q 


(4.5ii) 


(4.5iii) 


(4.5iv) 
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Here, we have a singular potential function (as well as a singular Green’s function). 

Now we shall touch upon the invariance of the potential equation (4.5ii I. In phase space 
the symmetry group is the group of canonical transformations, which preserve the canonical 
symplectic form. The group XO(3) is a subgroup of the canonical transformations. The 
required transformation is explicitly furnished by 


dpj + pjdq^ =dS{p, q ), 

S{p,q) := (c^' + Pi , 




Pj = 


dpj 

dS{-) 


=r jPk , 


(4.6i) 

(4.6ii) 

(4.6iii) 

(4.6iv) 

(4.6v) 


dqi 

n{q-,p) =n{q;p) , 

or, 

n (q; p) =n a) [q^ - c>) ; r)pk 

being the inverse transformation of and the parameters are in accordance with equa¬ 
tions ( |2.7i[ -iv)). Note that the particular transformation in ( |4.6iii| ) exactly corresponds to the 
group XO(3) of equation (2.7il. Therefore, the group XO(3) is also a symmetry group of the 
phase space potential Q. 


5 Discrete phase space: quantum mechanics and the potential function 

It is well known that in the one-dimensional idealized oscillator (reinstating the constants h 
and £ > 0), the energy eigenvalues are provided by 


2 [h'^ £2^ 


V’(Ar) “ ( ^ + 9 ] '^(N) > 


(5.1) 


Note that the operator -h has discrete eigenvalues in spite of P and Q both pos¬ 


sessing continuous spectra. The equation (5.1), interpreted phenomenologically, yields the 
discretization of the phase plane as 


ff £2 


1 

= n + -. 


(5.2) 


n G {0,1,2,-••} . 


The equation (5.2) represents denumerably infinite numbers of confocal ellipses in the back¬ 
ground of the continuous p — x phase plane. Discrete phase space in this section comprises 
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of three such equations 


1 

2 







(5.3) 


In other words, the particle’s accessible phase space is discrete in this case in such a way that 
the above combination of and pj is integer valued IJSI]. In much of the remaining section 
we again choose units so that the characteristic length £ = 1 and h= 1. 

Previously we had defined three partial difference operators in p.l| ), ( |3.2[ ), and p.3| ). Here 
we need to introduce two more partial difference operators in the following. The weighted 
mean difference operator is 


Af/(n^n^n^) := ^ 


\/+ !/(.., + 1 ,..) — — 1 ,..) 


and the weighted average difference operator is 

1 


Aj/(n\ n^) := ^ + 1,..) + - 1,..) 


(5.4) 


(5.5) 


Moreover /(n^, n^, n^) := 0 for any < 0. The physical significance of the operators (5.4 1 
and (|5.5|) is that they provide an exact discrete representation of quantum mechanics via 


PjxJj := — iAj n^) 


PjQ’^ - Q’^Pj 


xp = — i6j 'ip{n'^ ,rr). 


(5.6i) 

(5.6ii) 

(5.6iii) 


Now, we denote a Hermite polynomial, Hn{k) and define a related polynomial as follows: 






TT 


for e {0,1, 2,...} , 

poo 

I dkj S^j , 




(5.7i) 


(5.7ii) 

(5.7iii) 


(Here the index j is not summed.) We can infer from (5.7iil iii) that (kj) is the normalized 
eigenfunction of the momentum operator Pj = —iAj. 

The discrete phase space potential equation, analogous to the equations ( |2.1| ), is chosen to 
be 


- S^^PjPkW = y^A*A*W{n^,n^, n^) = 0 . (5.8) 
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A Green’s function for the partial difference equation (5.8 1 , which is analogous to that in (2.2 1 
is provided by: 


G(n^, n^, n^; n^, n^) := 


i=i 


[{k,Y + {k,Y + 




, (5.9ii) 

G(0, 0,0; 0,0,0) =2. (5.9iii) 

In analogy with ( 2.3i| ), we present the following potential from (5.9i I in this representation: 

IU(o)(n^, n^, n^) = G{n^,n‘^, n^; 0,0,0) 


(•jjn^+n^+n^ 


t 3/2 2("-^+"^+"^)/2(ni! n^! [(fci)^ + (fc2)^+(fc3)^ 


Iff 


g-(fef+fel+fe§) 



.2=1 


dkidk2dk^. (5.10) 


Restricting our analysis along the third axis for the momenj^(0, 0, 2n) admits the closed form 
expression: 


IU(o)(0,0,2n) = 


(-1)- ffr 

7r3/22nVSd JJJ [{ki^ + (^2)2 + (k^n 

— OO 

2^1+1^! 

--> 0, 

(2n + l)v^(2n)! 

IU(o)(0,0,2n + 2)l 
IU(o)(0,0,2n) 

1^(0) (0,0,0) = 2. 


dkidk2dks , (5.Hi) 


(5.11ii) 


2 n + 2 
2n + 3 


2 n + 1 
2 n + 2 


< 1 


(5.11 hi) 
(5.11iv) 


Thus, {lH(o)(0,0, 2n)}^ is a positive-valued, monotone decreasing sequence. Now we shall 
explore the asymptotic behavior in the limit re —)• oo VR(o)(0,0, 2re). For that purpose, we 


express (5.1 lii i as 


FU(o)(0,0,2re) = 


2n-ri + 1) 


(5.12) 


(2re + 1) y^r(2re + 1) 

The asymptotic expression for the gamma function is provided by 

r(re) = v/^(re) 2 e " [l + + 2 ^ 1 ? “ slsiiu “ 24883201? + ^ (l?)] ■ (5.13) 


^We will show later that the potential is rotationally invariant and so there is little loss of generality with this 
restriction and is done for mathematical convenience. 
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Therefore, equations (5.IIii), (5.12i, and (5.I3I imply that 


IU(o)(0,0,2n) = 


TT 


1/4 


f 1 

[1 + -! 
n 

n 


1 ^ 

4 + 


[(i+2k)='q 


1 + 12(n+l) + ^ 


({n+l)^) 


1 + 


+ 0 


24(2n+l) V (2n+l)^ J _ 
TT^/^ r 'I 

lim IU,'o)(0,0, 2n) =—0 • \/e • 1 • 1 > = 0 . 
n-^oo ' V® I J 

We plot fU(0, 0, 2n) in figurefor values of n near the origin. 


(5.14i) 


(5.14ii) 


2 . 0 | 

1 , 8 - 

1 . 6 - 

1.4- 

1 . 2 - 

1 . 0 - 

0 . 8 - 

0 . 6 - 

0.4- 

0 


Discrete phase-spaee potential W(0,0,2n) 


❖ 



o 


6 


<> 


o 

8 


❖ 


10 


Figure 1: The discrete phase space potential IU(0,0, 2n) for £ = 0.1. Note that this potential is finite 
everywhere and has a value of 2 at the origin. 


Now we want to compare the continuous phase space potential Q^,Q^',Pi,P 2 ,P 3 ) to 

the discrete phase space potential VU(o)(0) 0; 2n). Note that Qo(q^,q^, q^;pi,P 2 ,P 3 ) is a com¬ 
plex potential, and hence we really wish to compare the real, monotone VU(o)(0) 0) 2n) with 

the quantity y^|Qo(0) 0) 9 ^;Pi)P 2 ,P 3 )P- This comparison is also equivalent to comparing 
VU(o)(0) 0) 2n) to the potential V(o)(0,0, x^) = (l/47r)l/|x^| noting that the pj dependence 

drops out in the magnitude |Qo(0, 0, q^;pi,P 2 ,P 3 )l^- Hence, y^|Qo(0,0, q3;pi,p2,P3)|^ and 
V(o) (0,0, x^) differ only by a numerical factor of 1/\/2. 

From ( |5.3| ) we can use the approximation 

= (l/2f)(x^f-l/2 + 0(£^) (5.15) 


for extremely small values of £. Specifically £ is assumed small in comparison fo h, perhaps 
even as small as fhe Planck lengfh. Some caution musf be applied in fhe comparison as 
fhe £ and h are nof of fhe same dimension and hence fhe approximafion (5.151 only really 
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becomes valid for large values of x and provided the momentum is not large. In this limit the 
discrete phase space potential function essentially goes over to a configuration space potential 
function, making the comparison meaningful. Hence we limit our analysis away from the 
origin. (We have already shown in (S.llivI and figure |T] that the potential W(Q'^{'n}-, ti?, n^), 
without the approximations, is finite at the origin.) This comparison is done in figure|^for the 
various potentials. 


Comparison of potentials away from the origin 



Figure 2: A comparison of the potentials Vq (dashed), |f7(o) | (black) and lU(o) (diamonds, for ^ = 0.1) 
away from the origin. The discrete phase space potential is actually finite at the origin (see main text 
and previous plot) whereas the other two potentials diverge at the origin. 


The most general transformation of the scalar field W (n^ under the six-parameter 

(proper) continuous group XO+(3) is given 


W{n^ ,n‘^ ,rf) =exp 


3 4 


_y_U 

AjAf - AfcAf + Af A,- - Af A^ 




^kj = _ ^jk ^ 


lU(n^, n^, 
(5.16i) 

(5.16ii) 


The translation subgroup is parameterized by the continuous parameters and rotations are 
encoded in the continuous parameters The discrete phase space potential equation re¬ 
mains invariant since the generators and 


A,A 


# 

k 


^ 41 a ^ 

AkAf + A*A, 


Afh 

3 K 


of the group XO+(3) commute with the Casimir in¬ 


variant operator S^^A'^Af. Therefore the potential is also invariant under the continuous 
group XO{2>). 

Finally, we shall investigate the variational derivation of the discrete phase space potential 


^Here, 0''"(3) denotes the proper subgroup (rotation group) of the orthogonal group 0(3). 
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(5.8 1 . Consider a class of action sum 


ATI 7V| ATS 

'^[/] •= ^ ^ ^ 2 / 1 , 2 / 2 , • 


(5.17) 


The variational principle states: 


lim 
£—^0 


J[f + eh]-J[f] \ ^ 


The equation (5.18|) implies the following Euler-Lagrange equations: 


dL{- 


dy 


I- 


-A* 


an---) 

dVi 


= 0 


(5.18) 


(5.19) 


and boundary terms. The derivation of ( 5.18| ) and the explicit boundary terms will be shown 
in the Appendix. In the case of the discrete phase space potential equation, the appropriate 
Lagrangian is furnished by 


•= VU(n)AfIU(n) 


(5.20) 


(Recall that it is the momentum operator which yields the potential equation Laplacian and 

n O 

hence AJ (but not Aj) appears here.) 


6 Concluding remarks 

In this paper we have studied potential theory in discrete phase space . The Green’s functions 
are non-singular as are their associated potentials. Unlike the usual discrete theory which 
arises as an approximation of continuous systems, these Green’s functions have been shown 
here to be rotationally covariant. As well, the operators giving rise to these Green’s functions 
provide an exact representation of quantum mechanics. This provides a novel representation 
for quantum mechanics. It is known that singular Green’s functions, for example, lead to 
divergences in S'-matrix elements in relativistic quantum field theory. With the relativistic 
version of the representation presented here, the corresponding 5-matrix elements have been 
previously shown to be finite ifTOll . 


A APPENDIX - Variational principle in discrete phase space 


We will consider the two-dimensional lattice function /(n^, n^) for the sake of simplicity as 
generalizations to higher dimensions is straight forward. We use capital Roman indices to 
take values from {1,2}. The action functional is defined by 

Nl Nl 

Jif]~ Y. Y ^(y^^y^^y^\=fir.),yA=A*fir^)■ 

n^=Nl n2=Af2 


Subsequenfly, we denofe 
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The variational principle states that 


£—^0 


/ 


^2 

^l=A^l n2=Ar2 


^ h{n^,n‘^) + ^ - - Aj/i(n\n^) . (A.2) 

^2/ |(ni,n2) OyA |(nl,n2) I 


^2', 


We now use the following relation |8| 

,n‘^)A*h{n^ ,n‘^) + h{n^ ,n^)A^(f){v} ^v?) 


= Aa 


I -!,■) + h{r?,v?)(t){., - 1,.)] | , (A.3) 


where has been defined in (3.11. (In the above equation the summation convention has 
been suspended.) 

Next we utilize the discrete Gauss’ theorem ll^ for a discrete vector field 


ATI Nl 


n^=Nl n2=Ar2 


n 


(A.4) 


7V| 




[j\N^ + l,n^)-j\Nl,n^)]+ [f {n\ + 1) - jHn\ Nf)] 


^2=Nf 




Using ( A.3| ) and ( |A.4| ) in ( |A.1[ ), and utilizing the discrete Dubois-Reymond lemma |[24l . 
we derive the following: 


The Euler-Lagrange equations: 


dL{- 


dy 


-A 


# 


|(nl,n2) 


dL{- 


dVi 


= 0 


(A.5i) 


3 J 


for <n^ <N^ + 1 


The boundary terms: 

Nl 

E 


" ' + 1 


n2=Ar2 


dL{- 


/i(A^2 4“ 4“ 


dL{- 


h{N2,n^'' 


dyi |(Ari,n2) ’ dyi 


|(iVi+l,n2) 



+ 


2 

dL{- 




9yi \{Nl,n'2) 


-A 


#dL{---) 


Nl dL{- ■ ■) 


^y |(Afi\n2) dyA |(Arl,n2) V 2 5?/i |(iVi_l,n2) 



+ 1 


E 

ii=Ari 


dL{- 


dy2 |(ni,Af2) 


/l(?T. , N 2 + 1) + 


dL{---) 

dy2 |(ni,Af2+l) 


h{Nl,n^) 

/i(n\ A|) 
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Nl 


+ 


dLj---) 
dy2 |(ni,Ar2) 

-A 


/i(n^, Af — 1) 


dL{---) Nf dLj---) 

dy |(ni,Ar2) dyA \{n^,Nf) ^ 2 dy2 |(ni,Ar2-l) 


h{n^,Nf)'j = 0 . 

(A.5ii) 

Note that the discrete domain D 2 for the validity of the Euler-Lagrange equations is given by 
D 2 := {(n\ n^) : + 1, < r? < + l} . (A.6) 

However, the variationally admissible boundary of the discrete domain is provided by 
d*D 2 = {{n\n^) : + 1, - 1, 

U {(nS n^) : A^^ < < Aa^ + 1, - 1, A^^} . (A.7) 

The domain D 2 and the boundary d^D 2 are depicted in figure 


(A/, N^+1 ) 




• - 

(A/, A/) 


(A/, A/+7 ) 

-• 


(A/-7, A/) 



(A/+7,aA 


(A/, A/) 

•- 


(A/, Nj) 

-• 


• (A/+7,aQ 


(A/, A/-7 ) 


(A/, A/-7) 


Legend: 


boundary 
domain edge 

boundary and domain edge 


Figure 3: The domain and boundary for the discrete variational principle in a two dimensional lattice 
space. 
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